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Abstract. Interfacial fluctuation effects occuring at wedge and cone filling transitions 
are investigated and shown to exhibit very different characteristics. For both 
geometries we show how the conditions for observing critical (continuous) filling are 
much less restrictive than for critical wetting, which is known to require fine tuning of 
the Hamaker constants. Wedge filling is critical if the wetting binding potential does 
not exhibit a local maximum, whilst conic filling is critical if the integrated strength of 
the potential is attractive. This latter scenario is particularly encouraging for future 
experimental studies. 

Using mean-field and effective Hamiltonian approaches, which allow for breather- 
mode fiuctuations which translate the interface up and down the sides of the confining 
geometry, we are able to completely classify the possible critical behaviour (for purely 
thermal disorder). For the three dimensional wedge, the interfacial fiuctuations are 
very strong and characterised by a universal roughness critical exponent i/^ = 1/4 
independent of the range of the forces. For the physical dimensions d = 2 and d = 3, 
we show that the influence of the cone geometry on the fluctuations at critical filling 
is to mimic the analogous interfacial behaviour occuring at critical wetting in the 
strong-fluctuation regime. In particular, for d = 3 and for quite arbitary choices 
of intermolecular potential, the filling height and roughness show the same critical 
properties as those predicted for three dimensional critical wetting with short-ranged 
forces in the large wetting parameter {lo > 2) regime. 
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1. Introduction 

Following the seminal studies of Cahn and Ebner and Saam two fundamental 
questions concerning the nature of the wetting transition have drawn an enormous 
amount of attention 0, ^: Firstly, is the transition first-order or continuous? In 
particular, how does the order depend on the range and balance of the competing 
intermolecular forces? Secondly, for continuous (critical) wetting transitions, what role 
do the interfacial fluctuations of the unbinding fluid interface play in determining the 
critical exponents characterising the diverging lengthscales? 

The answers to both these questions reveal a host of reasons why the laboratory 
observation of fluctuation effects at critical wetting is beset with difficulties. In brief, 
the well developed fluctuation theory of wetting concludes P, |] that three dimensional 
critical wetting requires a fine tuning of the effective Hamaker constants appearing in the 
binding potential. This sensitivity means that critical wetting is a rather rare phenomena 
of which only very few examples are known for fluid- fluid interfaces 0, ^. Secondly, 
for systems with long-ranged forces, the critical exponents at critical wetting are mean- 
fleld-like (and depend on the range of the forces) and the interfacial roughness induced 
by fluctuations is very small. Even for short-ranged forces, where renormalization group 
(RG) theory based on effective interfacial models famously predicts dramatic fluctuation 
induced effects ^ , no appreciable deviations from simple mean-fleld like behaviour have 
been observed in either Ising model simulation |^ or recent experimental studies ||^. At 
the very least, this indicates that the asymptotic critical regime for critical wetting 
is much smaller than initially expected ||^. The conclusion inevitably forced upon us 
appears to be that fluctuation effects at wetting transitions are of negligible practical 
importance. 

The purpose of the present article is to compare and contrast these features of 
critical wetting transitions with the properties of continuous (critical) flUing transitions 
occuring in both wedge and cone shaped geometries. The thermodynamics of the fllling 
transition has been discussed independently by a number of authors |jlO|, |ll], [12| but only 



recently have the analogues of the above questions been addressed |]T3|, |T5|, [T^, [17 
Here we present full details of a novel interfacial theory of wedge fllling and also extend 
this approach to study fluctuation effects occuring at the fllling transition of a fluid 
adsorbed in a cone shaped geometry. These turn out to be completely different to 
those predicted to occur in the wedge. In particular we wish to draw attention to 
two striking features of conic fllling which have not been reported before. First, that 
the conditions for critical conic fllling, even more than conditions for critical wedge 
fllling, do not rely on any flne tuning of the binding potential structure and should be 
observable in cone geometries made from substrates that exhibit even strong flrst-order 
wetting. Secondly we point out what appears to be an unexpected connection between 
the geometry-induced aspects of critical conic fllling and the theory of the universal 



strong-fluctuation (SFL) regime of critical wetting 0, [T^. Speciflcally, from analysis 



of effective interfacial Hamiltonian models, we show that for the physically relevant 
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dimensions d = 2 and d = 3, and also for almost all choices of the intermolecular 
interactions, the influence of the conic geometry on the asymptotic critical behaviour and 
the detailed scaling structure of the interfacial height probability distribution function 
(PDF) as one approaches a critical cone filling transition is to mimic the way in which 
interfacial fluctuations determine the analogous exponents and the PDF for critical 
wetting transitions for planar wall-fluid (or fluid-fluid) interfaces belonging to the SFL 
scaling regime. This remarkable connection between conic filling and wetting and has 
already been noted in d = 2 from the explicit results of exact transfer matrix and 
replica trick studies of two-dimensional filling transitions (where, of course, the cone 
is analogous to a wedge) in both pure and impure bulk systems [0, Extending 
this theory our analysis shows that the fluctuation properties of three-dimensional conic 
filling are quite distinct from those occuring at three-dimensional wedge filling |jl6| and 
are instead related to some of the most dramatic predictions of the theory of wetting for 
planar wall-fluid interfaces. Thus in d = 3 the divergence of the equilibrium filling layer 
thickness and roughness at cone filling are universal and the same as the renormalization 
group predictions of Brezin, Halperin and Leibler for critical wetting with short- 
ranged forces in the large wetting parameter {u > 2) regime. 

Our article is arranged as follows. In the next section we briefly recall the well 
known fluctuation theory predictions for critical wetting detailing the critical behaviour 
predicted to occur in the SFL regime in d = 2 and d = 3. In §3 we present the 
main methods and results of the fluctuation theory of three dimensional wedge filling 
emphasising both the conditions for critical filling and the predictions for fluctuation 



dominated behaviour. A brief summary of these results has been given before |T^ 
We finish this section by turning to the special case of filling in d = 2. The role of 
a breather-like soft mode which increases the width of the interface as it translates 
up and down the wedge is shown to be directly responsible for the scaling behaviour 
of the PDF. In §4 we discuss three dimensional cone filling, at mean-field level and 
beyond, predicting universal geometry-dominated critical behaviour characterised by 
an appreciable interfacial roughness and critical exponents identical to those predicted 
for the planar wetting SFL regime. We argue that this equivalence of the PDF's at cone 
filling and SFL wetting is special to d = 2 and d = 3 only. We finish our article with a 
summary of our main results. 



2. Fluctuation regimes at critical wetting. 

Consider the preferential adsorption of a thick liquid film at the interface between 
a planar wall and a bulk vapour phase at chemical potential fi and a sub-critical 
temperature T. We suppose that at bulk two phase coexistence, corresponding to 
chemical potential fj, = /igat' wall- vapour interface is completely wet by the liquid 
phase for temperatures T > T^. Here we are interested in the values of the critical 
exponents describing the continuous divergence of the mean wetting film thickness 
Ztt ~ t'^^% roughness C,± ~ t'^'^^ and transverse correlation length ~ t'~'^^\ as 
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t' = (Ttt — T)/Tt, — > 0. The subscript vr is used to emphasise that these quantities 
are pertinent to planar wetting transitions. We denote the reduced temperature t' 
to avoid confusion with the filhng reduced temperature (q.v.). The standard theory 
of interfacial fluctuation effects at the wetting transition is based on analysis of the 
effective interfacial Hamiltonian 0, ^ 

H[l] = J rfx||(V/f + (2.1) 

where / is a collective co-ordinate describing the local height of the interface above 
the wall, S is the surface stiffness (tension) of the unbinding liquid-vapour interface 
and W{1) is the usual binding potential which accounts for the direct influence of the 
molecular forces. Note that throughout this article we shall suppose that the factor 
/?(= l/ksT) is absorbed into all the effective Hamiltonians, stiffnesses and potentials 
etc. Mean-field (MF) studies indicate that for critical wetting (with long-ranged forces) 
W{1) is necessarily of the form 

iy(/) = -^ + ^ + ...;/>0 (2.2) 

where a oc (T*^^ — T) and 6 > are effective Hamaker constants. Here the exponents p 
and q are specific to the particular ranges of the intermolecular forces. In d = 3 and for 
van der Waals forces the appropriate values are p = 2 and q = 3. For fixed dimension 
d < 3 the scaling theory of Lipowsky and Fisher |jl8| (which is supported by transfer- 
matrix and approximate non-linear renormalization group (RG) studies) predicts that 
the critical behaviour falls into three generic fluctuation classes labelled mean-field (MF), 
weak-fluctuation (WFL) and strong-fluctuation (SFL) scaling regimes. These regimes 
arise due to the interplay of the direct binding potential with the entropic repulsion 
originating from the interfacial wandering. We briefly recall the scaling argument here 
since the same ideas generalise to filling transitions in wedge geometries. First recall that 
for interfaces, the roughness and the transverse correlation length are related via 
the wandering exponent ^± ~ We write Ctt keeping the sub-script to avoid confusion 
with different wandering exponents that emerge for filling in wedges and cones. For 
thermal forces and c? < 3 it is well known that 

C = ^ (2.3) 



with finite for d > 3 and ~ vln^|| for the marginal case. Thus we are led rather 
naturally [0 to deduce that the bending energy term in ( p.l| ) and also the entropic 
effect of interface collisions with the wall gives rise to a fluctuation contribution to the 
effective potential of the form l~'^ with r = 2{d — l)/{3 — d) for d < 3. The competition 
between these three terms gives rise to the three fluctuation regimes mentioned above. 
In the MF and WFL regimes the leading term in 772 is still dominant and the wetting 
transition temperature is unchanged by fluctuations. The critical exponents are non- 
universal and are given by [|l8l 

/3r=^, u;^^=-^ ^r=c^r (2-4) 

q-p " ^{q-p) " 
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/3r=^, u^^ = -L±i- ^r=c^r (2-5) 

T — p " 2{t — p) " 

In the SFL, corresponding to p < r, and representing the critical behaviour in 
systems with short-ranged forces, the transition temperature is lowered by fluctuations 
and the critical exponents are universal. In d = 2 the critical exponents and the 
interfacial PDF, Pn{l), are very well known (see for example ^ and also references 
in 0, 3) and given by 

/5f^ = 1, = 1, z/|f = 2 and P^{1) = -^e"!^ (2.6) 

Notice that in both the SFL and WFL the interfacial fluctuations are very strong 
and ~ ^± for d < 3. The unusual bifurcation mechanism controlling the flxed point 
behaviour as — * 3 means that critical wetting in three-dimensions systems is rather 
different in the SFL regime, though for systems with long-ranged forces MF theory 
remains valid. Fluctuation effects are only predicted to occur for systems with short- 
ranged forces and are traditionally and most simply described by a bare binding potential 
of the form g |, 0, § 

(/) = -ae""' + fe-^'^' + ...;/> (2.7) 

together with the hard-wall restriction W{1) = oo for / < 0. Here n denotes the inverse 
bulk correlation length of the adsorbed bulk phase whilst the Hamaker constants a and 
b have the same interpretation as for long-ranged forces discussed above. The mean-fleld 
predictions for short-ranged critical wetting are therefore /3s = (In) and z/y = 1. The 
fluctuation theory of critical wetting with short-ranged forces was famously investigated 
by Brezin, Halperin and Leiblier using a linearised RG scheme and shows that the 
critical behaviour is strongly non-universal and falls into three regimes (labelled I - III) 
depending on the value of the dimensionless wetting parameter u where 



. - ^ (2.8) 

The three regimes correspond to the ranges u < 1/2, 1/2 < u < 2 and uj > 2 
respectively. Usually attention focuses on the value of the transverse correlation length 
critical exponent which shows the strongest non-universal critical behaviour: 



where the last case corresponds to an essential singularity. In all three regimes the 
roughness ~ y/T^ so that interfacial wandering has a much more pronounced influence 
on the equilibrium distribution of matter in the wetting fllm compared to systems with 
long-ranged forces. In regimes (I) and (II) the wetting temperature is unchanged from its 
mean-fleld value and the fllm thickness still grows logarithmically but with a different 
amplitude to mean-fleld theory. The most dramatic effect of interfacial wandering is 
manifest in the regime (III) which bears all the hallmarks of the SFL regime. Here the 
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wetting temperature is depressed by the fluctuations, and the critical exponents are 
given by 

(III) /3, = 1, z/x = ^ (2.10) 
with an equilibrium PDF which is a simple Gaussian. 

3. Fluctuation Regimes at Critical Wedge Filling. 

3.1. Three-dimensional wedge filling: phenomenology. 

Consider a three-dimensional wedge formed by the junction of two walls at angles 
±a to the horizontal with the height of the wall above the plane is described by a 
height function Zyj{x,y) = a\x\. Thus the wedge has a 'V shaped cross-section in the 
x-direction with the wedge bottom oriented along the y axis. Following our earlier 
discussion of wetting at planar walls we suppose that the wedge is in contact with a 
bulk vapour phase at saturation chemical potential n = fi^at- General thermodynamic 
considerations |jlO|, indicate that the wedge is completely filled by liquid for 



temperatures greater then the filling temperature Tp, which is specified by the elegant 
equality 

en{TF)=a (3.1) 

where 9Tr(T) denotes the temperature dependent contact angle of the liquid drop at 
the planar wall-vapour interface. Thus the influence of the wedge geometry is to lower 
the temperature at which the liquid-vapour interface unbinds from the wedge bottom 
compared to the planar wall. The filling transition occurring as t = {Tp — T)/Tp ^ 
may be first-order or critical corresponding to the discontinuous or continuous divergence 
of the equilibrium interfacial height (/q) = {l{x = 0, y)) measured from the bottom of the 
wedge. Critical filling is also characterised [0 by the divergence of distinct correlation 
lengths describing fluctuations in the interfacial height along {^y) and across {^x) the 
wedge and also by the roughness of the unbinding interface. Note that the cross- 
wedge correlation length essentially measures the width of the flat, filled region of the 
wedge. The flatness of the filled region simply reflects the absence of any macroscopic 
curvature of the meniscus at two phase coexistence. Recalling our earlier notation we 
define the wedge wetting critical exponents via 

(/o)~t-^", a-t-'^", (3.2) 

where we have adopted a suitable super-script, W (for wedge), for some of the exponents 
to avoid confusion with those pertinent to critical wetting and also for cone filling which 
we will consider later. 

Wedge filling can also be studied using a simple modification of the effective 
interfacial Hamiltonian and for rather open wedges corresponding to small a (where 



tana ~ a) the appropriate model is |13, 14, 15, 16 



H[l] = j dx!^^{Vlf + W{l-a\x\)'^ (3.3) 
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where l{x, y) denotes the height of the interface measured above the plane. As shown by 
Rejmer et al this model can be justified from analysis of a more general drum-head 
model and with a binding potential interaction depending on the local normal distance 
of the interface to the wall rather than the simple vertical height used in ( p.3| ). The 
simplifying approximations implicit in the above model do not effect any of the universal 
physics associated with the filling phase transition. 



3.2. Mean- field Theory 

At MF level the equilibrium interfacial profile follows from the solution to the Euler- 
Lagrange equation found from the minimization of (|3.3| ) 

ti=W'{l-a\x\) (3.4) 

which is solved subject to the boundary conditions that /(O) = dl/dx\Q = and 
l{x) ± ax as X ±oo. This equation can be integrated once and leads to 

simple condition for the mid-point height [0, |T^: 

^ = W{{lo))-W{l^) (3.5) 

Note that in the limit (/q) —>■ oo, the RHS —>■ ^6^/2 which follows from the 
Youngs equation in the present small angle limit. The above equation has an elegant 
graphical interpretation [|1^] from which one can determine the order of the filling 



phase transition occurring as T ^ Tp at fixed a or equivalently as a at fixed 

temperature. Thus the MF condition for first-order /critical wedge filling is that, at the 
filling temperature Tp, the binding potential does/does not possess a local maximum 
separating the minimum at Z = from the extremum at infinity. Thus wedges made 
from walls that exhibit critical wetting, show critical filling whilst wedges made from 
walls that exhibit first-order wetting show both first order and critical filling depending 
on whether the filling temperature is above or below the spinodal temperature Tg at 
which the first-order wetting binding potential develops a local maximum. Note that 
since the filling temperature can be lowered almost arbitarily by increasing the wedge 
angle a it follows that the wedge will typically show first-order filling if a is small 
and critical filling for larger values of a. The tricritical value of at which the order 
of the filling transition changes character is given by a* = 6tt{Ts) the precise value 
of which depends on how weak the first-order wetting transition is. Of course this 
argument presupposes that a spinodal temperature for the binding potential exists. 
This is certainly true for systems with short-ranged forces where it follows from the 
Landau theory of wetting. For example in the global phase diagrams of Nakanishi and 
Fisher the locus of spinodal temperatures (for different surface fields) corresponds 



to the analytic extension of the critical wetting curve. The prediction that wedges made 
from walls exhibiting first-order wetting show critical filling can certainly be tested in 
Ising model simulation studies. For experimental systems it may be that the mean-field 
condition for critical filling cannot be met because the spinodal temperature Tg is too 
low. However as we shall see, fluctuation effects lead to a second mechanism for critical 
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filling beyond the MF condition described above which may overcome this difficulty. 
Also, and perhaps of greater practical importance, we shall show that the mechanism 
for critical filling in a cone is rather different to that for a wedge and is much more easily 
fuimied. 

To continue, from ( |3.2| ) it immediately follows that at a critical filling transition the 
MF value of the order-parameter critical exponent is simply determined by the leading 
order decay of W{1). The exponent q plays no role and we find (3^ = 1/p which is totally 
different to the corresponding MF critical exponent (3s = ^/{q ~ p) for critical wetting. 
The MF values of the other critical exponents follow from analysis of the height-height 
correlation function (with y = y' — y) 

H{x,x';y) = {6l{x,y)6l{x',y')) (3.6) 

and its Fourier transform 

S{x,x';Q) = jdy e'^y H{x,x';y) (3.7) 

which exploits the translational invariance along the wedge. At MF level, S{x, x'; Q) 
follows from the solution to the Ornstein-Zernike equation 

dx"C{x, x"; Q)H{x", x'; Q) = 6{x - x') (3.8) 

where the direct-correlation function is defined by 

evaluated at the equilibrium height. Thus we arrive at the differential equation 

+ SQ2 + W"{l{x)-a\x\))S{x,x'; Q)=6{x-x') (3.10) 
To proceed it is best to consider the moment expansion 

oo 

S{x, x'; Q) = Y, Q^'^S^nix, x') (3.11) 

n=0 

and ffist concentrate on the zeroth moment Sq{x,x') which describes the position 
dependence of correlations in the fluctuations of the interfacial height across the wedge. 
This is found by using standard techniques to be 

H(xx') /•"li'^d^l'l^'l) dx 



So{x,x') = (^\i{x) \ - a) (|/(x')| - a) I . + 



2aW'{lo) S io (/(x 



a 



\2 



(3.12) 

where H{x) denotes the Heaviside step function {H{x) = l;x > 0,H{x) = 0; x < 0). 
Thus the structure of correlations across the wedge reflects the properties of the 
equilibrium height proffie l{x) found from the solution to the Euler-Lagrange equation 
( p.4|) , which is particularly simple. In essence, the interface is fiat such that l{x) ~ {Iq), 
for X < {lo)/a, while for x > {lo) /a the height decays exponentially quickly to its planar 
value Importantly, the lengthscale which controls the exponential decay is the wetting 
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Figure 1. A depiction of the coherent fluctuations of the interface along the wedge. 
The diagram shows the relevant lengths associated with the derivation of the effective 
one-dimensional Hamiltonian that describes such fluctuations 



correlation length ^y, which remains microscopic at the filling transition temperature. 
One consequence of this is that we can anticipate that the lateral width of the filled 
region of the wedge is trivally identified as ~ 2(/o)/« so that = ■ ^^^^ 
we conclude that So{x,x') is essentially position independent provided \x\, \x'\ < ^^f^ 
implying that the whole filled portion of the wedge fiuctuates coherently. Moreover that 
there is negligible correlation between height fiuctuations inside the wedge with those 
outside and simlarly between positions on either side of the filled region. From (|3.12| ) 
one can iteratively construct the position dependence of the higher moments using the 
integral equation 

S2n+2ix,x') = -T. J So{x,x")S2nix",x')dx" (3.13) 

which follows from the Ornstein-Zernike equation above. From this it follows that in the 
asymptotic critical (scaling) regime S2n{x,x') is position independent for \x\, \x'\ < ^^/S 
with 

^2.+2(x,x') = — ^ ( -M^Y (3.14) 



2W{{lo)) \W'{{lo)), 

and is negligible otherwise. This behaviour is consistent with the simple Lorentzian 
structure factor 

So{x,x') 
1 + Q^l 



^(^'^';Q) = ?^; \xlW\<U2 (3.15) 



where S'o(0,0) = a/2W'{{lo)) and identifies the MF correlation length along the wedge 
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Substituting for Iq we find the desired MF critical exponent 1/p + 1/2. The 

important observation here is that Uy ^ z/^ so that the fluctuations are highly anisotropic 
and dominated by the modes along the wedge. Indeed as one approaches the bulk phase 
boundary becomes arbitrarily larger than implying that the fluctuations become 
pseudo-one dimensional. The last lengthscale to calculate within the MF analysis is the 
midpoint roughness defined by = {[l{0, y) — (/q)]^) which can be identified from the 
relation 

ei = i/(0,0;0) (3.17) 

where the RHS can be obtained from the Fourier inversion of 5(0, 0; Q). This leads to 
the roughness relation for the d = ?> wedge 

i^-ik-X (3-18) 



which is one of the central results of the MF analysis. Substituting for the divergence 
of ^y and lo we find that all dependence on the binding potential exponent p cancels 
leading to a remarkable universal divergence 

~ t~~^ (3.19) 

valid for all ranges of intermolecular potential. In summary, the MF values of the 
lenghscale critical exponents at critical filling are 

Next we turn to the fluctuation theory of filling beyond MF approximation. 
3.3. Fluctuation Theory: breather modes and their excitations 

The MF analysis presented above is only valid if the fluctuations are small in some 
sense. Since the divergence of even within MF, is algebraic we can simply use the 



'contact condition' [|T8l that cannot be greater than (/q) to determine the Ginzburg 
criterion for the breakdown of MF theory. Thus we can assert that the above critical 
exponents are valid provided 

|^~r(^-^)<i (3.21) 

which in the limit t — implies that p must be less then 4. For p > 4, fluctuation effects 
dominate and we can anticipate (analogous to the WFL and SFL regime behaviour of 
critical wetting) that the roughness is comparable with the interfacial height (Zq). 
Unfortunately it is very difficult to develop a full fluctuation theory of filling transitions 
based on the effective interfacial Hamiltonian ( p.3|) since, as is clear from the above 
remarks, the upper critical dimensions for wetting and filling are quite different. The 
situation is rather similar to that encountered in the theory of wetting where one would 
ideally like to develop a fluctuation theory based on analysis of, say, a semi-infinite 
Ising or Landau-Ginzburg Wilson (LGW) model. Such a approach would afford a 



Wedge filling, cone filling and the strong-fluctuation regime 



11 



description of both bulk (and surface) critical behaviour and also of wetting. Since this 
has not proved possible one instead focuses on effective models constructed to capture 
the essential physics of the phase transition in question. Of course, all Hamiltonians are 
effective Hamiltonians but some are more effective than others. The LGW model based 
on a local magnetization order-parameter can (in principle) describe bulk and surface 
criticality in addition to wetting and filling. In turn the effective interfacial Hamiltonian 
( p. 3D based on a collective co-ordinate l{x, y) is a theory of wetting and filling. The task 
here is to develop a simpler model than (|3.3|) which just affords a description of filling. 

The theory of wedge filling which we propose is constructed from the full interfacial 
model under the assumption that the only fluctuations that determine the asymptotic 
critical behaviour arise from the pseudo-one dimensional long-wavelength undulations 
in the local height of the filled region. We refer to fluctuations in the local height as 
breather modes since these increase/decrease the local width of the filled region (and 
the volume of liquid beneath it). Thus we seek to derive a model based on a collective 
co-ordinate lo{y) which only describes these modes (see Fig 3.1) and which is valid 
for sufficiently small wave- vectors Q ^ ^/ix- Such a model can be easily constructed 
by generalising the method introduced by Fisher and Jin p, Q for the derivation 
of wetting interfacial models from the LGW Hamiltonian. Assuming that for a given 
constrained configuration of the mid-point height l^iy) that all other fluctuations are 
small (relative to the soft mode), it follows that the desired filling Hamiltonian can be 
identified by 

eM-HwU) = Tr{eM-Hm (3.22) 

where the Tr denotes a partial trace with respect to the generalised crossing-criterion 
(GCC) constraint that l{x = 0,y) = lo{y). Since this trace is over non-critical 
fluctuations one can employ a saddle-point idenifiction 

HwMy)]=H[l^{x,y-lo{-))] (3.23) 

where /^(x, /o(')) denotes the profile that minimises ( |3.3|) subject to the GCC. As in 
the Fisher- Jin theory of wetting, the desired non-planar constrained configurations can 
be constructed perturbatively in terms of constrained profiles that are translationally 
invariant along the wedge (analogous to the planar constrained magnetization profiles) 
which we simply write as l^{x] Iq) to avoid a proliferation of indices. These translationally 
invariant profiles satisfy the Euler-Lagrange ( p.4|) subject to the GCC /^(x = 0; Iq) = Iq. 
The resulting wedge filling model has the form 

HwM = ldy[^{^y + ^^(^o)| (3.24) 

where Vw{lo) denotes the wedge binding potential and cr(/o) is a position dependent 
wedge line stiffness. Formally these quantities are determined by the relations 

Vw{h)= I dxl^i'^^^j^] +W{l\x-k)-a\x\)\ (3.25) 
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(up to unimportant additive constants defined such that Vw{Itt) = 0) and 

,(,).Sr.J*i^V (3.26) 



Sio 

The planar constrained profiles /^(x; /q) have the same simple near- fiat structure as the 
equilibrium MF profile discussed earlier. This means that the evaluation of ( p.25 , 3.26| ) 
is particularly simple so that for first-order and critical wedge filling we can approximate 



HwM = dy{^{-^) + Vwilo) > (3.27) 




The qualitative and quantitative form of the filling potential depends on the order of 
the MF phase transition. For critical filling it has the expansion 

W(/o) = ^^^'-""'^^^ + apll-P + ... (3.28) 
a 

and of course fluctuations are restricted to the region Iq > analogous to the hard- 
wall condition for wetting. Here the effective Hamaker constant ap = cl/{p — !)« and 
note that the coefficient of the linear term is simply proportional to t, the relevant 
reduced temperature. The critical filling potential shows a single minimum (located 
precisely at the MF value of Iq) which continuously moves out to infinity as t 0. 
The one-dimensional filling Hamiltonian can be studied using a number of different 
techniques including transfer-matrix theory |T5|, |16|, |2^ and approximate RG methods 
|23| . However there is no need to go to these lengths to extract the fluctuation regimes 
and the values of the critical exponents since they follow from a simple extension of 
the Lipowsky-Fisher scaling theory considered earlier. To see this we suppose that the 
influence of fluctuations at critical filling follows from analysis of an effective potential 

V^jj{lo) = Vwilo) + Vfl{lo) (3.29) 

which accounts for the entropic repulsion of the unbinding interface from the wedge 
bottom. Again the fluctuation term can be estimated from both the form of the bending 
energy contribution to ( 3.27 ) and also the number of collisions the interface has with 
the wedge bottom. This we estimate to of order Zq where denotes the thermal 
wandering exponent for wedge filling (in d = 3) relating the roughness and pertinent 
correlation length ~ The value (p = 1/3 follows from considering the invariance 
of the free part of the filling Hamiltonian under the scale transformation y ^ y/h and 
lo ^ b-^^lo. Thus we arrive at an effective filling potential of the form 

V ^ff{l,) = Vwilo) + CI f (3.30) 

which has a single minimum at the equilibrium mid-point height. In this way we find two 
different fluctuation regimes for filling depending on whether the direct intermolecular 
potential or entropic contribution determines the next-to-leading order correction to the 
linear term in ( |3.28D . The leading linear term is always relevant (in the RG sense) and 
reflects the influence of the wedge geometry on the statistical/thermodynamic properties 
of filling. This is rather similar to the theory of fluctuation effects at complete wetting 
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Figure 2. Schematic illustration of the wedge filling potential on approaching a first- 
order and continuous filling transition. 



transitions In this way we predict there exists a filhng MF regime (FMF) for 

p < 4, with exponents given by ( p.2CI| ) and a filhng fiuctuation dominated regime (FFL) 
for p > 4 with 

/?f = 7, ^y = \^ ^±=7; (3.31) 



4' - 4' - 4' 
As a final remark we turn to the effective one-dimensional Hamiltonian description of 
first-order filling. This is based on the model ( p.27|) but with a different potential VV(/o) 
the form of which is shown in Fig. 3.3. In fact the structure of Vw^o) for this case 
can be qualitatively deduced by noticing that the extrema of VV(^o) must correspond 



to the possible solutions of the MF equation ( |3.5|) . For Q.„ > a, Vwih) has only a 
single minimum which remains finite at T = Tp. However for 6^^ < a a potential 
barrier appears separating the local minimum and the infinitely deep potential well at 
infinity. Note that exactly at the first-order filling transition temperature the minimum 
of Vw{lo) has a finite well-depth relative to Vw{oo). This is quite unlike the behaviour 
of W{1) at first-order wetting transitions. These remarks indicate even if the mean- 
field condition for critical filling is not fulfilled, the infiuence of the one- dimensional 
breather mode fiuctuations may still drive the transition critical by causing the interface 
to tunnel away from the local minimum of the potential. Whether this happens or 
not depends on the range and depth of the minimum of Vw{lo) but will certainly be 
most important for systems with short-ranged forces (such as Ising models). Note this 
tunnelling mechanism does not alter the filling transition phase boundary (|3.1|) which 
simply refiects the relevance of the linear term in VV(^o)- One may also note that this 
tunnelling must definitely occur out of coexistence for temperatures T > Tp where it is 
responsible for the washing out of the pre-fiUing line. 

We mention only briefiy here that it is trivial to investigate critical wedge filling 
from off-bulk co-existence using the effective model by allowing for a term of order P 
which couples to the undersaturation in the filling potential. This clearly represents the 
contribution from the volume of undersaturated liquid in the filled region of the wedge. 
From this one can readily establish that the filling thickness, roughness and correlation 
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lengths all exhibit scaling with a gap exponent 6 = 5/4 in the FFL regime. Thus 
the divergence of the mean interface height along the critical wedge isotherm as one 
approaches saturation is characterised by a —1/5 power law. 

3.4- Wedge Filling in two dimensions 

It is a simple matter to extend the wedge filling fluctuation theory to arbitrary 
dimensions d and in particular d = 2 where one can compare the predictions with 



the exact transfer matrix solution of the full interfacial model ( |3.3| ). Consider a wedge 
in bulk d dimensions which is of 'V shape cross-section in the x direction but which 
is translationally invariant in the remaining d — 2 (surface) dimensions. Denoting the 
position of the interface above the wedge bottom by where x\\ denotes the {d — 2) 
dimensional vector parallel to the wedge orientation. The filling Hamiltonian trivially 
generalises to 

Hwik] = J dx^i 1^ (V/o)' + VV(/o)} (3.32) 

which clearly describes the {d — 2) dimensional thermal excitations of the breather 
interfacial mode. The scaling theory also generalises in a natural way. First we calculate 
the d dimensional wedge wandering exponent (w by considering the invariance of the 
free part of ( p.32[) under scale transformations. This determines (w as 



4- d 

Cw = ^ (3.33) 

relating the roughness and correlation length measured along the wedge (which we here 
denote ^y) via ^± ~ ^^'^ . Similarly the fluctuation contribution to the effective potential 
arising from the bending energy and entropy terms scales as 

Vf{Io) ~ ^o"^ = (3.34) 

where tw is the wedge fluctuation exponent analogous to r in the Lipowsky- Fisher 
theory of wetting. 

This heuristic scaling theory predicts that for general dimension d < 4, the 
criticality at filling falls into a mean-field (FMF) and fluctuation dominated (FFL) 
regime depending on whether p is less or greater than the marginal value p* = 
2{d — — d). The critical exponents are given by 

-,-4. -r~ + ^; FMF (3.35) 

p p 2 4 2p 



and notice that it is only in = 3 that the roughness exponent is universal independent 
of the range of the forces. These predictions for critical singularities at wedge filling are, 
in general very different to the predictions for critical wetting. However for the special 
case of (i = 2 there appears to be a connection between filling and critical wetting 
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discussed below. In d = 2 the breather fluctuation theory predicts that the marginal 
intermolecular range is p = 1 and that the values of the FFL regime critical exponents 
are 

= z/f = 1; d = 2 FFL (3.37) 

with Uy undefined. Note that these are pertinent to all physical examples of 
intermolecular forces. The first thing to observe here is that the numerical values 
of the filling critical exponents (3^ and uY in the d = 2 FFL are identical to the 
analogous critical wetting exponents in the {d = 2) SFL regime. As we shall now see 
this coincidence extends to the full form of the PDF. 

As well as predicting the critical exponents at filling the fluctuation theory also 
predicts the position dependence of the scaling form of the mid-point PDF, Pw{lo)- To 
see this note that the case of wedge filling in d = 2 is special because the breather mode 
is only subject to zero dimensional fluctuations since there is no direction along the 
wedge. As a result one can go beyond the above predictions for the values of the critical 
exponents and assert that in the asymptotic critical region, the mid-point PDF is given 
by 

P^{Iq) = Are-^^('«) (3.38) 

where here (and below) denotes a suitable normalization factor. Now from ( |3.28|) we 
have Vf(/o) ~ ^(^^ ~ a'^)lo/a ~ 2S(0^ — a)lQ so that in the critical regime the breather 
fluctuation theory leads to the scaling form of the PDF 

P^^(/o) = Ne-^^is.-^)io (339^ 

or equivalently 

Pwik) = Ae-T4> (3.40) 
{l-o) 

which is of course consistent with the values of the critical exponents quoted in ( |3.36| ). 
This prediction is remarkable for two reasons: (I) it is in exact agreement with the 
transfer matrix result for the PDF based on the {d = 2) interfacial model ( |3.3| ). For 
quite arbitrary choices of binding potential the mid-point interface height distribution 
function is given by 0, |l^ 

Pwik) = Ne'^''^'\¥^\lo)\' (3.41) 

where \E'o'^''(/) is the ground state eigenfunction appearing in the transfer matrix 
spectrum for (i = 2 wetting at a planar wall. This satisfies the well known Schrodinger 
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- ^§2'^o\i) + Wmi;'\l) = Eo^i^'^l) (3.42) 

where Eq is the corresponding groundstate eigenvalue which is related to the contact 
angle by Eq = — /2. From the properties of "^^^^ (l) it follows that the filling transition 
is indeed characterised by the scaling regimes and corresponding critical singularities 
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predicted by the general fluctuation theory. Moreover for p > 1, and far from the wall 
the wave function decays as 

^W(/)^e"^'-' (3.43) 



so that ( p.39| ) and ( |3.41| ) are identical. The equivalence of the scaling form of the 



PDF with the breather mode fluctuation theory and the full effective interfacial model 
completely supports our conjecture that it is the overall breather mode which translates 
the interface up and down the sides of the wedge, rather than the thermal wandering of 
the interface within the filled region that determines the criticality at filling. 

(II) The derived expression for the filling PDF in the FFL regime has an identical 
structure to the analogous interfacial height PDF, Pnil) for critical wetting in the SFL. 
Recall that this quantity is determined as Pjril) ~ |\I'o'^''(/)p with the transfer matrix 
approach so that we in the SFL we have the well known result, quoted earlier 

PAl) = f (3.44) 

which is identical to the result for filling. In fact in ci = 2 the connection between filling 
and critical wetting for purely thermal disorder extends to the case marginal case p = 1 
for filling where the structure of the FMF/FFL PDF has the same form as the PDF for 
critical wetting at a WFL/MF borderline for which the specific heat exponent vanishes 
15| . The PDF for this case is not simply a pure exponential function. 

The equivalence of the film thickness (and roughness) critical exponents and PDF's 
for two dimensional, fluctuation-dominated filling and wetting has also been established 
for the case of filling/wetting with random bond disorder described by the interfacial 



model 15 



H[l] = J ||(£)2 + V;(x,/(x)) + H^(/-«|x|)| (3.45) 

with disorder averages satisfying 

V{xj{x))=0 (3.46) 

and 



Vr{x, l{x))Vr{x', l'{x)) = A5{x - x')S{l{x) - l'{x')) (3.47) 

where A is a measure of the disorder. The study of filling transitions in this model allows 
us to test the possible connection between fluctuation dominated two dimensional filling 
and critical wetting in a system where the interfacial wandering exponent is different 
to the thermal case. It transpires that for systems with purely short-ranged forces 
the model ( |3.45|) can be solved exactly using the transfer matrix method by extending 



Kardar's Bethe Ansatz approach |^ used in his analysis of wetting. This again confirms 



the thermodynamic prediction for the location of the filling transition and identifies the 
universal filling critical exponents as 

= 2, uY = 2 (3.48) 
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which are identical to the analogous critical exponents (3s and ux, determined by Kardar 
for critical wetting with random bonds with short-ranged forces corresponding to the 
SFL regime 0, |2^, . The expression for the scaling form of the mid-point interfacial 
height PDF Pf{Io) for d = 2 random bond filling is highly non-trivial but again analysis 
shows that it is identical to the PDF for random bond wetting [jl5|. 



4. Conic Filling. 

4.I. Phenomenology. 

Consider a non-planar wall-fluid interface in the shape of an infinite cone which makes an 
angle a to the horizontal. Thus, for open wedges (small a) the height of the wall above 
the plane is described by a height function zc{r) = ar with r the displacement vector 
along the two-dimensional displacement vector along the horizontal plane with origin 
at the cone apex. As before the wall is supposed to be in contact with a vapour phase 
at two-phase coexistence with thermodynamics indicating that the cone is completely 
filled for T > Tp, where the conic filling temperature Tp satisfies the same condition 
^ff(Tp) = a as filling in the d = 3 wedge (^.11). We allow for the possibility of first-order 
and critical cone filling corresponding to the discontinuous or continuous divergence 
of the mid-point filling height (/q) as t = (Tp — T)/Tp 0. Also of interest is the 
roughness arising from the interfacial fluctuations of the unbinding liquid vapour 
interface. The radius of the near fiat, filled region of the cone, which also characterises 
height fluctuations in the central region is trivially related to the height via ~ (^o)/*^- 
As with wedge filling we anticipate that for radial positions r > the interface will 
simply take the microscopic value it would have on a planar wall. Thus we define two 
lengthscale critical exponents for critical cone filling by 

(/o)~t-^% (4.1) 

in an obvious notation. Here we seek to understand both the conditions under which 
critical cone filling can occur and the values of the critical exponents at MF level and 
beyond. As we shall show these are quite distinct from the situation for wedge filling. 

Following our earlier treatment of wedge filling we anticipate that interfacial 
fluctuations manifest themselves in two different ways at conic filling. Firstly there is 
the usual wandering of the interface within the filled region determined by the value of 
the standard wandering exponent Ctt which is of course marginal in c? = 3 (concentrating 
on thermal fluctuations). Such fluctuations would occur even if the edges of the filled 
region were fixed to a given height, /q say, and would contribute terms of order ^Iq and 
Vln /q to the interfacial roughness in ci = 2 and d = 3 respectively. However much more 
important than this is the breather mode that translates the entire flat region up and 
down the cone. The fundamental distinction between three-dimensional wedge and cone 
filling arises because the breather mode is thermally excited by quasi-one dimensional 
and zero- dimensional excitations respectively. In this respect three dimensional conic 
filling is the higher dimensional version of two-dimensional filling considered in the last 
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section. Note that the effective dimensionahty of the breather mode does not destroy 
the phase transition because the ffuctuating field diverges as the filhng temperature is 
approached. 



4-2. Mean- field theory. 

Denoting the local interfacial height above the plane by /(r), the natural starting point 
for the study of conic filling for open cones (small a) is the standard interfacial model 

H[l]= [ df\^ [Vlf + W{l-ar)\ (4.2) 



. 2 

where, again, r denotes the radial vector from the apex. To begin we chose binding 
potentials W{1) of the type (|2.2| ) corresponding to walls that exhibit critical wetting 
and set the exponent parameter q = p — 1. Numerical minimization of the model 
Hamiltonian for different ranges of the intermolecular potential shows that the filling 
transition is critical (for the present choice of critical wetting walls) and located precisely 
at = tt- However the numerics also indicates that the MF value of the critical 
exponent characterising the divergence of (/q) appears to fall into two distinct regimes. 
For p > 1 the critical behaviour is universal and dominated by the geometry such 
that f3^ = 1. On the other hand if p < 1 the critical exponent is consistent with the 
identification (3^ = 1/p, identical to that found within the MF theory of wedge filling. 
This intriguing MF critical behaviour follows from a very simple and highly accurate 
variational theory of conic filling. Noting the simplicity of the equilibrium interfacial 
profiles obtained in our numerical minimization we parameterise the possible height 
function by 

= I ^ ^ 7? (4.3) 

^ ^ 1 + ar, r > ^ ^ ^ 

and determine the mid-point height dependence of the conic filling potential defined by 

VSo) = H[l^]-H[l^] (4.4) 

Note that this variational approximation becomes increasingly accurate for as one 
approaces a critical filling transition since the profile tends to a macroscopically flat 
meniscus. Minimization of the conic filling potential leads to the following integral 
equation for the mid-point height 

E{a^-9l){lo~Q no 
2 

This equation has an simple graphical interpretation which recovers the numerical 
results quoted above and leads to an elegant criterion for the order of the filling 
transition. The filling transition is continuous provided 



W{l)dl (4.5) 



W{l)dl < (4.6) 
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corresponding to an integral measure of the overall attractive/repulsive nature of the 
binding potential. By allowing for pertubations to the simple variational ansatz it is 
easy to show that this is the exact MF condition for criticality within the effective 
interfacial model to order 0{a^). The condition for criticality of conic filling is much 
less severe than that for both critical wetting and critical wedge filling and implies that 
even substrates which exhibit rather strong first order-wetting will exhibit a continuous 
filling transition in a cone geometry. A comparison of the respective shapes of the 
binding potential need for criticality for three dimensional wetting, wedge filling and 
cone filling is shown in Fig. 4.1. Provided the condition for critical cone filling is met 
the MF critical exponent falls into two regimes: 

/3f = ^ for p < 1 and /3f = 1 for p > 1 (4.7) 

The universality of the critical behaviour in the geometry dominated critical regime 
is reminiscient of the way in which fluctuation effects lead to universal critical behaviour 
at wetting transitions. This connection is explored further below. 



4-3. Fluctuation theory: breather modes and the PDF. 

The effective zero- dimensionality of the breather mode at conic filling suggests that we 
should analyse fluctuation effects in the same manner as we did for filling of a wedge/cone 
in d = 2. Recall that this simple approach gave the precise scaling form of the mid-point 
interfacial height PDF (and hence the critical exponents) at filling as that found from 
the exact transfer-matrix theory. Accordingly we propose that the critical behaviour at 
conic filling is described by a PDF for the mid-point interfacial height is of the form 

Pc{lo) = Are-^<^('«) (4.8) 

since the role of the cone potential in d = 3 is analogous to the wedge potential Vw{lo) 
in d = 2. It is a simple matter to show that V^(/o) has the form 

Veil) = " - crcl + act' + ■■■ (4-9) 
a 

consistent with the MF values of the critical exponents discussed above. The three 
terms in this expression represent the purely geometrical contributions arising form the 
area and circumference of the filled region and also the direct influence of the specific 
intermolecular forces. Both ac and Oc are positive effective Hamaker constants. For 
p > 1 the final term in ( [4.8| ) is irrelevant, geometry dominates and the PDF can be 
written 

Pcik) = Ne-^.''-'-''' (4.10) 
and identifies the universal values of the critical exponents as 

= 1 and z/^ = ^ for p > 1 (4.11) 
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Figure 3. Schematic illustration of the qualitive form of the binding potential W(V) 
necessary for a continuous wetting/filling transition. For critical wetting at a planar 
wall (top) the potential decays monotonically exactly at T^. For a wedge, the potential 
must show no local maximum at Tp. For a cone the integrated stengh of the potential 
must be attractive {A^ < A~'). 



On the other hand if the second term in (|4.9| ) is irrelevant and close to its most probable 
position the PDF is Gaussian with critical exponents given by 



1 J C 1 r 

and z/7 = - for p < 1 



(4.12) 



p - 2 

Equations ( ^.6| ) and ( [4.8[ - [4.12D are the main predictions of this section. Here we make 
two pertinent remarks: 

(I) the MF values of the filling height critical exponent (3f' is unchanged by 
fluctuations. However the thermal fluctuations of the breather mode do lead to an 
appreciable interfacial roughness. The predictions ( [4.11| ) are appropriate for both 
experimental systems with long-ranged dispersion forces and simulation studies of Ising- 



Wedge filling, cone filling and the strong-fluctuation regime 



21 



like systems with short-ranged forces. For discrete lattice based Ising studies, it is 
probably more convenient to study adsorption in an inverted pyramid shaped geometry. 
The critical exponents are the same as the cone for this case. 

(II) the scaling form of the PDF and the values of the critical exponents in the 
universal geometry dominated regime of cone filling are identical to the analogous critical 
singularities characteristic of the SFL regime of critical wetting short-ranged forces i.e. 
regime (III) corresponding to a; > 2. This is directly analogous to the observations 
made earlier for wedge/cone filling in d — 2. 

The above predictions point to a remarkable coincidence between the universality 
of conic filling in d — 2 and d = 3 and the SFL regime of critical wetting in the same 
dimension which is encapsulated in the equivalence of the film thickness and roughness 
critical exponents as well as the interfacial height PDF's. Thus in the important 
geometry dominated regimes of cone filling and for the physically relevant dimensions 
d = 2 and d = 3 one can identify 

Ps = = ^f"- and Pcik) = Pr\lo) (4.13) 

4.4- Universal conic filling and the SFL regime. 

The above predictions beg the question, why does the influence of the cone geometry 
on the (universal) unbinding at filling mimic the fiuctuations characteristic of critical 
wetting in the SFL? Usually, of course, if two different systems exhibit the same critical 
exponents one concludes, after identifying the appropriate order-parameter, that they 
belong to the same universahty class. We do not beheve that this is the correct 
explanation because in this case there does not appear to be an analogue of the wetting 
correlation length critical exponent i/y for filling. In other words the equivalence of the 
one-point PDF's does not extend to the two point function. In view of this we forward a 
simple plausibility argument suggesting that the equivalence of the PDF's is not simply 
coincidence but is special to the physical dimensions d = 2 and d = 3. We emphasise 
that the specific results of our model calculations for the values of the critical exponents 
stand quite apart from this discussion. 

Let us focus on the structure of the PDF for critical wetting at a planar wall 
in d — 2 and d — 3. The distribution function P7r(0 found by evaluating the 
(normalised) Boltzmann sum over all interfacial configurations that pass through the 
specified height I at an arbitarily chosen position f — (say), equivalent to the 
identification P^(Z) = Ne~^*^''^ where F* — —InZ^ and Z* is the partial partition 
function summing over this class of configurations. In turn the constrained frcc-cncrgy 
F*[l] is a functional of the constrained height profile corresponding to the (partial) 
average over all profiles that satisfy the appropriate condition l{x = 0) = /q. For 
large / ^ 1^^ the constrained profile determining the PDF naturally forms the shape 
of a triangle (in d = 2) and cone (in d = 3) characterised by the equilibrium contact 
angle 9t^. Thus the droplet shape of the constrained profile determining the asymptotics 
of PttII) near critical wetting is essentially identical to the equilibrium height profile 
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measured with respect to the local wall height near critical cone filling (since 9^^ ^ a). 
This observation hints at a possible connection between the PDF's at wetting and cone 
filling. One may also reason that if there is any connection between the PDF's of filling 
and wetting then it will only be valid when the direct effect of intermolecular interactions 
is unimportant since the non-universal influence of these forces manifests itself in very 
different ways. This is only true in the geometry dominated regime of cone filling and 
the SFL regime of critical wetting. 

On the other hand it is easy to see that any connection between the two PDF's 
is not valid for all dimensions. The generalisation of our simple fluctuation theory of 
cone flUing presented to arbitary dimension d is immediate. Using the same variational 
approach described above we arrive at a cone binding potential of the form (neglecting 
unimportant constants) 

Vc{l) = tl'^-' - f-"" + • • • (4.14) 

representing the leading order surface and line contributions. The influence of 
intermolecular forces is present through a term 0{l'^^^^^) which is higher-order for 
p > 1 independent of the dimensionality. Using this cone potential to evaluate the PDF 
arising from the (zcro-dcimcnsional) thermal fluctuations of the breather mode leads 
to the following prediction for the universal, geometry dominated cone flUing critical 
singularities in arbitrary dimension 2 < d < 4, 

/3f = 1 and = (4.15) 

which of course recovers our results for wedge/cone filling m. d = 2 and cone filling 
in (i = 3 as special cases. For > 4 the roughness remains finite. Only for the 
physically important dimensions d = 2, 3 to the numerical values of these critical 
exponents coincide with the predictions for SFL regime critical wetting. Nevertheless 
we find it remarkable that for these special dimensions not only the values of the critical 
exponents but also the detailed form of the PDF's are identical for fiUing and critical 
wetting. Finally we mention that the exponents quoted above allow us to identify the 
value of the wandering exponent for critical cone filling defined by ~ Noting 
that interfacial height is trivially related to the radial correlation length (width) of the 
filled region we are led to the identification 
4 - d 

Cc--^ (4.16) 

which shows the infiuence of the zero-dimensional breather mode on the distribution of 
matter. Obviuously in two dimensions Cc = Cvy = 1 since the cone and wedge geometries 
are identical. 

5. Conclusions. 

In this article we have presented details of a general fiuctuation theory of filling 
transitions at wedges and cones. The main results and conclusions of our study are 
summarised below: 
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(I) The dominant interfacial fluctuations at wedge and cone filling correspond to 
breather- like motion that translates the flat, flUed region up and down the sides of the 
conflning geometry. 

(II) For three dimensional wedge filling the thermal excitations of the breather 
mode are described by a novel one-dimensional effective Hamiltonian which predicts 
a universal value of 1/4 for the roughness critical exponent for all ranges of the 
intermolecular forces and also large scale interfacial roughness and universal critical 
singularities for p > 4. 

(III) For fiUing in d — 2, the simple picture of an effective zero-dimensional 
fiuctuating breather mode correctly captures the critical exponents and precise scaling 
behaviour of the mid-point height PDF as found from the full transfer-matrix solution 
to the full interfacial model (of wetting and filling). 

(IV) For two-dimensional filling with both thermal and random bond disorder the 
scaling of the mid-point height PDF is identical to the corresponding results for SFL 
regime critical wetting in this dimensionality. 

(V) Conic filling shows two different types of critical behaviour including a universal 
geometry dominated regime for p > 1 pertinent to all physical types of intermolecular 
force. Applying the simple breather-like fluctuation picture to this phase transition 
yields a scaling form for the PDF and associated interfacial height and roughness 
critical exponents identical to the SFL regime prediction for critical wetting in this 
dimensionality. We have argued that this equivalence of the PDF's is special to c? = 2 
and d = 3. 

(VI) Even at MF level the conditions under which three dimensional wedge and 
cone flUing are continuous are different to the flne tuning required for critical wetting 
and are sensitive to the particular geometry. For critical wedge wetting one requires that 
the binding potential shape at T — Tp shows no local maximum. On the other hand 
cone-flUing is critical provided the integrated weight of the potential remains negative 
(attractive). 

These conclusions, regarding both the conditions for criticality and the 
manifestations of interfacial fluctuations at wedge and cone filling, are in sharp contrast 
to the situation for critical wetting. Concentrating on the values of the thickness and 
roughness critical exponents as pertinent to realistic three-dimensional systems with 
non-retarded van der Waals forces {p = 2) we have 



indicating that the geometry has a profound infiuence on the roughness of the interface. 
This reflects the special role of the breather mode at flUing transitions. Consequently if 
critical wetting at a wall-fluid interface were ever seen experimentally, the observation of 



/?s = 1 and i^x = (In) for critical wetting 

1 1 

— - — 7 for critical wedge flUing 

(3^ — 1 and i'^ — ^ for critical cone flUing 



(5.3) 



(5.1) 



(5.2) 
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wedge and cone filling in this system would certainly provide a means of observing large 
scale interfacial fluctuations. It should certainly be possible to test the predictions 
of universal and large scale interfacial fluctuation related behaviour in Ising model 
systems. We have also emphasised that the geometry dramatically alters the conditions 
for criticality and certainly relaxes the requirement of fine tuning the values of the 
Hamaker constants necessary for critical wetting. The condition for inducing continuous 
filling transitions is easiest to fulfil for the cone geometry and is highly encouraging 
that the observation of critical filling in cones made from walls exhibiting first-order 
wetting is experimentally feasible. It may even be possible to further manipulate the 
conditions for criticality by considering filling in other geometries and by introducing 
chemical heterogeniety. Throughout this article we have focused exclusively on the filling 
transition occuring as T — Tp at bulk two-phase coexistence. Also of interest is the 
infiuence of the geometry on complete wetting adsorption isotherms corresponding to 
the continuous divergence of the filling thickness for temperatures T > Tp as /i fi^^^^ 



|28| , P9| . These transitions have also received much recent attention and show great 
sensitivity to the shape of the confining geometry. For cones and wedges in all dimensions 
the divergence of the filling height is super-universal and given by 

/o ~ ^ and a ~ ^ (5.4) 

where 5fi denotes the undersaturation. However unlike critical filling the breather mode 
does not change the nature of the roughness so that ~ It would be very 

interesting to explore the possible connection between conic filling and SFL critical 
wetting in two and three dimensions in more detail. Whilst for d = 2 the equivalence 
of the critical exponents and PDF's has been establised exactly using transfer-matrix 
methods it would be extremely instructive to understand this from a RG perspective. 
This would perhaps throw light on a more detailed theory of three dimensional cone 
filling based on the full interfacial model. Also our analysis of wedge and conic filling 
in d = 3 was limited to pure thermal disorder unlike our analysis in ci = 2. Studies of 
three dimensional wedge and cone filling in impure systems would be very illuminating. 
It may be that for impure systems, filling and SFL critical wetting only share common 
features in d = 2. 
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